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ABSTRACT 

Current sheets thinner than the ion inertial length are unstable to the tearing instability and will 
develop magnetic islands that grow due to magnetic reconnection. We investigate whether the growth 
of magnetic islands in a current sheet can continue indefinitely, or in the case of the heliosheath 
until reaching a neighboring current sheet, and at what rate the islands grow. We investigate the 
development and growth of magnetic islands using a particle-in-cell code, starting from particle noise. 
Performing a scaling of the growth of magnetic islands versus the system size, we find that the growth 
rate is independent of the system size up to the largest simulation we were able to complete. The 
islands are able to continue growing as long as they merge with each other and maintain a high aspect 
ratio. Otherwise there is not enough magnetic tension to sustain reconnection. When applied to the 
sectored magnetic fields in the heliosheath, we show that the islands can continue growing until they 
reach the sector width and do so in much less time than it takes for the islands to convect through 
the heliosheath. 



At low latitude the solar wind is divided by the he- 
liospheric current sheet into sectors of oppositely di- 
rected magnetic fields. The th ickness of th e current 
sheet, A, is around 10,000km (jSmithl l200lh at 1AU, 
and the separation between each sector, or the sec- 
tor width, is around 1AU. The sector width remains 
nearly constant, although increasingly variable, out to 
the termination shock (TS) at ~ 9 0AU, where the super- 
sonic solar wind ab ruptly slows (Beh annon et al.l 1989; 
iBurlaga et al.ll2007f ). For the essentially collisionless en- 
vironment of the solar wind A controls whether collision- 
less reconnection onsets. For A greater than the ion in- 
ertial scale, di = c/uj p ^ where c is the speed of light 
and Co> p i is the ion plasma frequency, reconnection via 
the collisionless tearing ins tability does not take place, 
while for A < di it does (jCassak et al.l 120051 : lYamadal 
120071 ). Just upstream of the TS where the ion density is 
~ 0.001cm -3 , di ~ 7200km, which remains smaller than 
A based on measurements at 1AU. Thus, the absence 
of significant reconnection of the sector field upstream of 
the TS is consistent with models, although a definitive 
study of A upstream of the TS remains to be carried out. 
A downstream of the TS is predicted to be 2500km based 
on measurements at 1AU and the shock compression, 
while di is 4200km based on ion density m easurements of 
about 0.003cm- 3 (jRichardson et al.ll2008l ). Thus, down- 
stream of the TS the current sheets should begin break- 
ing up into magnetic islands. 

The growth of islands in the finite plasma j3 (the 
ratio of the plasma pressure to the magnetic pres- 
sure) heliosheath has been investigat ed in parti cle-in- 
cell simulations dDrake et al.l 120101 : IQpher et all 120111 : 
iSchoeffler et al.ll2011f ). The predicted island length in 
the initial phase of reconnection is around 190^, much 
smaller than the sector width, 8900di, so it is uncertain 
whether islands will grow to the full sector width. Due 
to computational limitations it is not possible to real- 
istically simulate the disparate ion inertial and sector 
scales. In this work, we perform a scaling study of the 



growth of islands, in which the inter-current sheet sep- 
aration (equivalent to the sector width) varies, in order 
to understand what happens in the real system. Even 
our largest simulations have sector widths that are much 
smaller than in the actual sectored heliosheath. 

The current sheets separating the sectored regions be- 
gin to form islands after crossing the TS. As the islands 
grow, the current sheets are convected toward the he- 
liopause. The plasma flows outward at around 80km s _1 
and steadily decreases in speed for 20AU, at which point 
the r adial flow remains close to zero (|Krimigis et al.l 
|2011[ ). The important question which we seek to answer 
is whether the islands are able to expand to the sector 
width before the current sheet reaches the heliopause. If 
the islands expand to the sector width, there would be no 
more laminar field that can shield cosmic rays, and cos- 
mic rays could percolate through the system consisting 
solely of islands. In addition, the full expansion of the 
islands would imply that acceleration mechanisms due to 
the growth of islands may play a significant role. 

We simulate the growth of magnetic islands using a 
particle-in-cell code, p3d. The initial conditions con- 
sist of two op positely dire cted current sheets in Har- 
ris equilibrium (|Harrislll962f ) . with a superimposed back- 
ground density. The initial magnetic fields are in the x 
direction, which corresponds to the azimuthal direction 
in the heliosheath. The current flows in the z direction, 
which corresponds to the north south direction. The y 
direction corresponds to the radial direction of the he- 
liosheath. In the heliosheath the islands are predicted 
to be highly elongated due to the development of pres- 
sure anisotropy, and this elongation is dependent on both 
the io n-to-electron mass rati o and the electron temper- 
ature ([Schoeifler et al.l l2011f ) . Typical simulations use 
a reduced mass ratio in order to reduce computational 
expenses, which produces much shorter islands than ex- 
pected for the real system. We therefore use an enhanced 
temperature of the electrons in the background in or- 
der to form more realistic elongated islands. In the he- 
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Fig. 1. — Out-of-plane current density, j z , for t = (a) 60f2 ci , (b) 
80^, (c) 120n^ 1 , and (d) 160^ 1 For better contrast, all points 
with \j z \ > 1 are assigned the colors shown for 1 or —1. 

liosheath nearly all of the pressure comes from the ions 
due to the pick-up ion population. In our simulation 
half of the pressure is in the ions and half is in the elec- 
trons. The /? of the plasma in the asymptotic field is 3, 
which corresponds to a temperature in the heliosheath 
of 650,000K, based on a density, n = 0.003cm -3 , and a 
magnetic field, B = 0.15nT. The ratio of the proton to 
electron mass in this simulation is 25. We have run simu- 
lations that imply the reconnection rate is not sensitive to 
the value of the mass ratio, so our prediction should scale 
to the realistic mass ratio. The background temperature 
is 15 times the Harris sheet temperature of 0.2brriiC 2 Al 
where ra$ is the ion mass and ca is the Alfven speed, 
based on the asymptotic magnetic field, Bq. The ratio 
of the speed of light to the Alfven speed is 25. Each sim- 
ulation has a grid scale resolution of = A y = 0.05di 
and a time resolution of dt = 0.0040^, where Q C i is 
the ion cyclotron frequency. The half thickness of the 
current sheet is set to wo — 0.5di so that collisionless 
reconnection can begin from particle noise. We simu- 
late a two-dimensional system. In three-dimensional sys- 
tems islands form at different z locations and grow along 
z, eventually stagnating likely due to interaction s with 
other islands (Sha y et alJ2QQ3tlSchreier et al.ll20l(l . The 
significance of this effect in the heliosheath is unknown. 

The time evolution of the largest of these simula- 
tions, with dimensions of 409. 6di x 102.4^, can be seen 
in Figure [TJ By t = 601^ 1 the current sheet has 
broken into elon gated magnetic islands as predicted in 
iSchoeffler etHI (|2011f ) (Figure Ufa)). The length of the 
islands is smaller than the separation between the two 
current sheets, so it is expected that they could not grow 
to the sector width since circular islands do not have ten- 
sion to drive reconnection. However, as can be seen in 
the subsequent times (Figure QJb)-(d)), the islands on 
a given current sheet begin to merge. Merging length- 
ens the islands which enables further growth until they 
approach the neighboring current sheet. 

The islands found in the heliosheath, which are pre- 
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Fig. 2. — Width of the largest island in the lower of the two 
current sheets vs. time for simulations with dimensions of 102 Adi x 
25.6<2i (green), 204.8^x51.2^ (red), and 409.6^x102.4^ (black). 
The plus signs denote the point where the island width reaches 44% 
of the system size and represent the time when the island begins 
to interact with the other current sheet. The blue dashed line is a 
line of best fit between these three points. The slope of this curve 
is the rate of growth at which the island expands, 0.12cAt>- 

dieted to be much shorter than the current sheet sepa- 
ration, should in principle grow to the sector width as 
long as the islands are able to continue merging. The 
aspect ratio of our simulation box is 4, which was suf- 
ficient for islands to continue merging up to the time 
that they begin interacting with the neighboring current 
sheet. We also carried out a simulation with aspect ratio 
2 where the islands stopped merging when their length 
approached the system size. The islands saturated and 
stopped growing when there was no longer enough mag- 
netic tension to maintain reconnection. The size of the 
heliosheath in the azimuthal direction greatly exceeds the 
sector width. Islands should therefore grow to the sector 
width. The next question is how long it takes for these 
islands to reach the sector width. In order to establish 
this time we carry out a scaling of the time required to 
grow to the sector width for system sizes 102.4^ x 25. 6<^, 
204.8di x 51.2di, and 409.6d* x 102.4^. 

To determine the size of the island we find the mini- 
mum of the flux function , ip where the magnetic field, 
B = z x (x, y) + B z (x, y) z, along the center of the 
lower initial current sheet at a particular time. This min- 
imum corresponds to the most developed x-point. The 
upper current sheet has x-points at the maxima. The 
line of constant flux that crosses an x-point is known as 
the separatrix. The distance between the maximum and 
minimum ^-locations of the separatrix is defined as the 
island width, w. Figure [2] shows the island width ver- 
sus time for each of the simulations. We use the time 
when the island reaches 44% of the size of the box as a 
measure for when the island reaches the neighboring cur- 
rent sheet. At 50% the island begins to be affected by 
the presence of the neighboring current sheet. The best- 
fit line connecting the times when the islands reach the 
neighboring current sheet fits very well with the island 
width versus time for all the simulations. 

The results of Figure [2] suggest that the islands grow at 
a nearly constant rate that is independent of the system 
size. Keeping in mind that this number is based on only 
three data points, using the slope of the best fit curve 
we can obtain an estimate for the rate of growth of the 
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Fig. 3. — Reconnected flux in the lower of the two current sheets 
vs. time for simulations with dimensions of 102. 4c^ x 25.6^ (green), 
204.8d; x 51.2di (red), and 409.6^ x 102.4^ (black). The plus signs 
denote the point where the island width reaches 44% of the sys- 
tem size and represent the time when the island begins to interact 
with the other current sheet. The blue dashed line is a line of 
best fit between these three points. The slope of this curve is the 
reconnection rate 0.079Boc^. 

island of around 0.12c^6, where CAb is the Alfven speed 
based on the background density, which is distinct from 
the normalization of the code, c^, which uses the peak 
density of the Harris equilibrium. If we extrapolate the 
growth rate to a very large system we can predict a time 
for the islands to reach the sector width. Based on a 
magnetic field, B, of 0.15nT and density, n, of 0.003cm -3 
the Alfven speed just downstream of the TS is 60km s _1 . 
Using this speed for c^6, and the sector width, W = 
0.25AU, we obtain a growth time, t g , of about 60 days, 
much less than the plasma convection time across the 
heliosheath. 
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Assuming the radial velocity of the solar wind inside the 
heliosheath is 70km s _1 during island growth, this time 
corresponds to a distance of 2.5AU past the TS. 



In addition to the rate of change in the island width 
being nearly constant, we find that the reconnection rate 
is also independent of system size. In Figure [3] we look 
at a similar plot to Figure [2j but of the reconnected 
flux. To find the reconnected flux we take the difference 
between the maximum and minimum of the flux func- 
tion along the center of the initial current sheet. Based 
on the slope of the best fit curve, the reconnection rate 
was 0.079. Previous scalings of reconnection rate versus 
system size for conditions with lower plasma /?, relevant 
to the Earth's magnetosphere and the 1AU solar wind , 
have shown comparable rates (|Shav et al.l I1999L [2007). 
However it is i mportant to note t hat t he reconnection 
rates shown in iShav et al.l (|1999l I2007D are associated 
with a steady state reconnection, as opposed to the rate 
shown here which includes the complicated dynamics of 
the merging process as well. 

In the range of sizes simulated the rate of flux recon- 
nection and the rate of island growth are nearly con- 
stant once reconnection begins. The growth rates are 
independent of the system size. It is reasonable to con- 
clude that in a larger system these trends would con- 
tinue. The merging of magnetic islands allows the islands 
to maintain a high aspect ratio, which maintains the 
magnetic tension necessary to drive reconnection. The 
steady reconnection rate allows for a constant rate of is- 
land growth, resulting in islands with a width that scales 
like the current sheet separation. These islands would 
be fully grown long before reaching the heliopause. The 
growth of these islands in the heliosheath is vital for the 
generation of anomalou s cosmic rays (AC Rs) by Fermi 
acceleration in islands (Drak eet al.l l2Q10h . Since these 
islands are expected to be present in the sectored region, 
and the flux o f ACRs is greatly re duced outside of the sec- 
tored region (Op her et al.ll2011| ). both observations and 
models suggest that the sectored heliosheath has broken 
into magnetic islands. 
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